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We present a unified description of the relativistic nNN and jnNN systems where the strong 
interactions are described non-perturbatively by four-dimensionai integrai equations. Our formu- 
lation obeys two and three-body unitarity and is strictly gauge invariant with photons coupled to 
all possible places in the strong interaction model. At the same time our formulation is free from 
the overcounting problems plaguing four-dimensional descriptions of nNN and ~/nNN systems. The 
accurate nature of our description is achieved through the use of the recently introduced gauging of 
equations method. 

I. INTRODUCTION 

Recently we have introduced a method for incorporating an external electromagnetic field into any model of hadrons 
whose strong interactions are described through the solution of integral equations The method involves the 

gauging of the integral equations themselves, and results in electromagnetic amplitudes where an external photon is 
coupled to every part of every strong interaction diagram in the model. Current conservation is therefore implemented 
in the theoretically correct fashion. Initially we applied our gauging procedure to the covariant three-nucleon problem 
whose strong interactions are described by standard four-dimensional three-body integral equations More 
recently we used the same method to gauge the three-dimensional spectator equation for a system of three-nucleons 
11- 

Here we apply our gauging procedure to the more complicated case of the covariant nNN system whose four- 
dimensional integral equations have only recently been derived |^,^| . In this way we obtain gauge invariant expressions 
for all possible electromagnetic processes of the nNN system, e.g. pion photoproduction 7c? — > nNN, pion electro- 
production ed —* e'nd, and because pion absorption is taken into account explicitly, we also obtain gauge invariant 
expressions for processes like deuteron photodisintegration -fd — » NN and Bremsstrahlung NN — ► <yNN that are 
valid even at energies above pion production threshold. Included in the electromagnetic processes described by our 
model is the especially interesting case of elastic electron-deuteron scattering. Here the deuteron is described as a 
bound state of our full irNN model; thus, after gauging, we obtain a rich description of the electromagnetic form 
factors of the deuteron with all possible meson exchange currents taken into account. 

II. FOUR-DIMENSIONAL nNN EQUATIONS 

Since the early 1960's many attempts have been made to formulate few-body equations for nNN -like systems using 
relativistic quantum field theory |7|-|i"^] . Yet all these attempts have had theoretical inconsistencies with perturbation 
diagrams being either overcounted or undercounted. The first consistent nNN equations were derived by us only 
recently Q and have the feature of containing explicit subtraction terms in the kernels of the equations that eliminate 
all overcounting. At the same time the undercounting problem is eliminated by the retention of certain three-body 
forces. The same nNN equations were later derived by Phillips and Afnan Q using a different method. In this section 
we would simply like to restate these equations but in a form that is particularly convenient for gauging. 

A. Distinguishable Nucleon Case 

It is easy to rearrange the four-dimensional nNN equations of Ref. 5 into a convenient form similar to the one used 
by Afnan and Blankleider |L3| in a three-dimensional formulation of the nNN system. For the distinguishable nucleon 
case we obtain 
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FIG. 1. Illustration of quantities making up V^. (a) V^ Ed , (b) B, and (c) F 2 = fcd^ 1 . 



T d = V d + V d G?T d (1) 
where T d , V d , and Qf are 4x4 matrices given by 

rrd _ ( T NN T$ \ vd _ f V NN T A r d _ ( D A , , 



T d T d ) ' " ~ ^ T d G^l J' yt ~ \ G oW °G 

Here we use the conventions of Ref. 13 and refer to the two nucleons as particles 1 and 2, and the pion as particle 3; 
also A = 1, 2 labels the channel where nucleon A and the pion form a two-particle subsystem with the other nucleon 
being a spectator, and A = 3 labels the channel where the two nucleons form a subsystem with the pion being a 
spectator. The superscript d is used to denote the distinguishable nucleon case on those symbols that will later be 
used without the superscript for the case of indistinguishable nucleons. Eq. ([l]) is a symbolic equation representing a 
Bethe-Salpeter integral equation to be solved for T d . T d consists of transition amplitudes T^ N , Txn, 2V m , and Ta m 
(A and \x are the spectator-subsystem channel labels) the last three being elements of the matrices T$, T$, and T d , 
respectively. The physical amplitudes for NN—>NN, NN—>Trd, nd—^NN, and Trd—^ird are then given by 

Xnn = t nn ! X dN = i'dTsN ; X Nd = T N3 ^p d ; X dd = -ipdT^-iipd, (3) 

respectively, where ipd is the deuteron wave function in the presence of a spectator pion. The kernel V d specified in 
Eq. (0) consists of the following elements: 

Vnn = vSn™ + B d G B d - A (4) 

where V®^ Ed is the nucleon-nucleon one pion exchange potential illustrated in Fig. |l|(a), B d = B + PBP where B 
is the diagram illustrated in Fig. |l](b) and P is the nucleon exchange operator, and A is a subtraction term that 
eliminates overcounting. T d is a 3 x 1 matrix whose A'th row element is given by 

2 

Ti = Y J h i F i -B d . (5) 

»=i 

Here S\i = 1 — 8\i and Fi — fidj 1 (i,j = l or 2 with i ^ j) where is the iVj — > nNi vertex function and dj is the 
Feynman propagator of nucleon j. F2 is illustrated in Fig. |l](c). Note how B d here plays the role of a subtraction 
term. T d is the 1x3 matrix that is the time reversed version of T d (similarly for other "barred" quantities) , Go is the 
ttNN propagator, and X is the matrix whose (A, /i)'th element is Sx.^. Finally the propagator term Qf is a diagonal 
matrix consisting of the NN propagator Do, and the 3x3 diagonal matrix w° whose diagonal elements are tid^ 1 , 
t%di , and tfd^ 1 , with tx being the two-body t matrix in channel A (for A = 1 or 2, tx is defined to be the nN t 
matrix with the nucleon pole term removed). 

The subtraction term A is defined with the help of Fig. || as 

A = + W d N + W NN + X + Y d (6) 

where W d N = W nN + PW^ N P and Y d = Y + PYP. 
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FIG. 2. Terms making up the subtraction term A. (a) W-mr, (b) W^n, (c) Wnn, (d) X, and (e) Y. The dark circles 

represent the following two-body amplitudes: (a) full tttv t-matrix, (b) one-nucleon irreducible nN t-matrix, and (c) full NN 
t-matrix minus the NN one-pion-exchange potential. 
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B. Indistinguishable Nucleon Case 



In approaches based on second quantisation in quantum mechanics it is usual to obtain the scattering equations 
for identical particles by explicitly symmetrising the equations of the distinguishable particle case. This procedure is 
not justified in the framework of relativistic quantum field theory. Nevertheless, as we have already derived the irNN 
equations taking into account identical particle symmetry right from the beginning ||, we can formally deduce how 
the above distinguishable nucleon irNN equations need to be modified in order to get the indistinguishable nucleon 
case. With this in mind, we introduce identical nucleon transition amplitudes defined in terms of distinguishable 
nucleon transition amplitudes as 

Tnn = Tfj N A _ Tan = Tin A TdN = T 3 nA 

Tna = T N i — T N 2P T N d = T N3 A Taa = Tu — Ti 2 P (7) 

TdA = I31 - T 32 P TAd — Ti 3 A T d d = T 33 A 

where A — 1 — P is the antisymmetrising operator. As Tni = PTn 2 P and Tn 3 = PTn 3 P, we can alternatively 
write Tna — ATni, TNd — ATn 3 , and TdA = AT 3 \. Thus in the transition to the indistinguishable particle case, 
the original 16 transition amplitudes for distinguishable particles have been reduced to 9 antisymmetrised transition 
amplitudes. By taking residues of the irNN Green function for identical nucleons at two-body subsystem poles, one 
obtains the following expressions for the physical amplitudes: 

Xnn — Tnn ', XdN — tpdTdN ', XNd = TNd^d ', Xdd — ipdTddipd- (8) 

Using Eqs. (0) in Eq. (0) it is easy to show that one again obtains a Bethe-Salpeter equation 

T = V + VG t T (9) 

but where now T, V, and Qt are 3x3 matrices given by 

Tnn T N a T N d \ ( Vnn^ AT\ F 3 A \ 

T = ( T AN T A a T Ad ; V = TfA -G^P G^A ; (10) 

T dN T dA T dd J V Fi A G o' A / 
\D \ 

Qt=\ Gohd^Go (11) 
lG t 3 d 3 x G J 

where t 3 = t 3 A is the t matrix for two identical nucleons in the presence of a spectator pion. 
From Eq. (0) it follows that 

VnnA = V^ E - W* - - W&f - X R - Y LR + l -B LR G B LR (12) 



where VS^ E = V^ Ed A and where we have used a superscript L to indicate that A is acting on the left, and 

7 R = W L W R * T = Wb„ 

7T7T V " 7T7TI V 'NN V 'NN 



superscript R to indicate that A is acting on the right. It is easy to see that W R ^ = W^, W R N — Wh[ N and 



X R = X L . 

We note that the 7riYAT equations for identical nucleons derived in Ref. 6, unlike Eqs. (0), ( [Toj ) and (|TT|) , contain no 
antisymmctrisation operators A and are therefore not equivalent to our irNN equations. 



III. ttNN electromagnetic transition currents 



A. Gauging the irNN Equations 

In this section we shall derive expressions for the various electromagnetic transition currents of the irNN system. 
To do this we utilise the recently introduced gauging of equations method . As the gauging procedure is identical 
for the distinguishable and indistinguishable particle cases, we restrict our attention to the ttNN system where the 
nucleons are treated as indistinguishable. 

Direct gauging of Eq. (||) gives 

t^ = v^ + v^g t T + vg?T + vg t r^ (13) 
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which can easily be solved for T M giving 



T" = (1 + rOt)V(l + GtT) + TQ?T. 



(14) 



7~ M is a matrix of gauged transition amplitudes T NN , T NA , T Nd , etc. To obtain the physical electromagnetic transition 
currents of the nNN system where photons are attached everywhere it is not sufficient to just gauge the physical nNN 
amplitudes of Eq. (|^). Although this would indeed attach photons everywhere inside the strong interaction diagrams, 
it would miss those contributions to the physical electromagnetic transition currents where the photons are attached 
to the external (initial and final state) pions and nucleons. In order to also include these external leg contributions it 
is useful to attach the corresponding propagators to the A-amplitudes of Eq. (£3): 



Xnn = DqXnnDo 



XdN — d^XdNDo 

Xdd = dixXdd d n - 



(15) 
(16) 



The physical electromagnetic transition currents are then obtained by gauging Eqs. (^|) and (^) and "chopping off" 
external legs: 



JNN — u 
JNd ~ ^0 A Nd u tt 



- d ^ 1X dN D o 
d^X^d- 1 . 



3dd 



Using Eqs. (@) we obtain that 



3mr = D o 1d o t nn + T NN D%D 1 + T NN 
JdN = 4>d D T dN + d^ 1 (j)dG^TdN + 4>dD T% N 

+ 4> d D Q T dN D%DQ l 
JNd = T Nd D Q ^ + TNdG^dd' 1 + T£ d D (j>d 

Jdd = (t>d D o T ddD (j)d + d~ X (j)dD%TddD 4>d 
+ 4>dD T^ d D (j)d + (frdDoTddD^dd' 1 
+ (^dDoTddDo^d 



(17) 
(18) 



(19) 
(20) 
(21) 

(22) 



where <f>d is the deuteron bound state vertex function defined by the relation ij) d = d%d2<pd- Note that (f) d consists of 
contributions where the photon is attached everywhere inside the deuteron bound state, and is determined by gauging 
the two-nucleon bound state equation for <fi d [§],3 ■ 

That the above equations are gauge invariant is evident from the fact that we have formally attached the photon 
to all possible places in the strong interaction model. The gauge invariance of our equations also follows from a strict 
mathematical proof; however, as this proof is essentially identical to the one given for the NNN system ||], we shall 
not repeat it here. 



B. Alternative Form of the nNN Equations 



Although the preceding discussion solves the problem of gauging the nNN system in a straightforward way, the 
expression obtained to calculate the gauged transition amplitudes, Eq. (|l4|), is not the most convenient for numerical 
calculations. The disadvantage of Eq. ( [L4| ) is that it utilises a Green function Q t which contains two-body t matrices. 
The presence of these t matrices makes the evaluation of T M unnecessarily complicated. We therefore present an 
alternative form of the nNN equations which uses a "free" Green function which contains no two-body interactions 
and which leads to simpler expressions for the nNN electromagnetic transition currents. 

The nNN equations Eqs. (^|), ( |l0| ) and (pi]), can be written in the form 



f Tnn 










=( 



Vnn 
T 











Tnn Tn 
T n T 



(23) 



where Vnn = V NN A is given by Eq. (p"2|), and 



4 



T d A \ [ (F 2 - B d )A 
T$A J ~ [ (Fx + F 2 - B d )A 

T = {AT( Ajf) = [A(F 2 - B d ) A{F X + F 2 - B a 



-P A 
A 



M^ 1 




|M 3 1 



(24) 
(25) 
(26) 



With the view of gauging Green function versions of irNN transition amplitudes, we introduce the Green function 
matrix appropriate to Eq. (123): 



G = 



Gnn Gn 
Gn G 



AD 




Gnn Gna GncI 
= | Gan Gaa GAd 
GdN GdA Gdd 



D 
Go 



Tnn Tn 
T n T 



D 
Go 



(27) 



The inhomogeneous term is chosen so that Gnn corresponds exactly to the full Green function for AW scattering. 
Then it can be shown that Q satisfies the equation 



G = 



where A and A are given by 



A 



AD„ 
LG 



- \B d ) A 



1 



\V NN - |ALG A |A 
|A t 



A=[A (F - \B d ^ 



\B d 



The essential feature of Eq. (B8h is that it is written in terms of an effective "free" Green function matrix 



Go = 



AD 
LG 



(28) 



(29) 



(30) 



which does not involve two-body interactions. For this reason Eq. (28) is ideal for the purposes of gauging. 



C. Gauging the Alternative Form of the nNN Equations 

In terms of the elements of G, the Green function versions of the physical amplitudes (defined in Eq. (|l6|)) are given 

by 

X N n — Gnn X d N = ^dGdN X N d = GNd^d X d d — ^dGdd^d- (31) 
After gauging, these equations give 

X NN = G NN X dN = <j)^GdN + 4>dG dN (32) 

X Nd = G Nd (t>d + G N d4>d x dd = 4>dGdd(t>d + 4>dG^ d 4> d + 4>dG d d4>d- (33) 



The irNN electromagnetic transition currents j^a are then determined by Eqs. (|17| ) and Eqs. (|ii 

To determine the quantities G^ in Eqs. (^) and (^3[) we need to derive the expression for by gauging Eq. (|2~g|). 
Defining 

Vt = / \V N n - jAiGoA |A\ ^ 

Eq. can be written as 

G = Go + GoVtG- (35) 

Gauging this equation and solving for G^ gives 
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FIG. 3. Graphical representation of (a) and (b) which enter into the expression for the gauged Green function 
matrix of Eq. (|4(i|). Constant factors and effects of operator A and matrix Q. have been ignored in this illustration. 



g» = {i + gv t )g%(i + v t g) + gv?g. (36) 

To simplify this equation we cannot use Eq. ( |35|) to write 1 + Vt5 = g^ X g and 1 + gVt = gg^ 1 since L is singular so 
that the inverse g^ 1 does not exist. Instead we use the fact that L — LQL, where 

n = ! 2 O (37) 
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which allows us to write 



where 



QS = GoM"g„ (38) 



_ ( i^sr'rSDf ) (39) 

Using Eq. ( |38| ) in Eq. ( |36| ) gives us a compact expression for 

0" = g (A-P + VH g. (40) 

It is easy to see that the term corresponds to photon coupling in the impulse approximation while Vf corresponds 
to the interaction currents and consists of the elements = (Vnn — ALGoA)^, A M , A M and V*. It is important to 
note that the diagonal elements of matrix i M are both of the form 

(Ud- 1 )" = ttdf + U (dj 1 )" = t^J 1 - ( 41 ) 

where — d~ l d^d~ l is the electromagnetic vertex function of the nucleon and the last equality follows from the fact 
that (dj 1 dj) tl = 0. Thus the diagonal elements of t M involve new subtraction terms UY^ whose origin does not lie in 
the subtraction terms of the strong interaction ttNN equations, but rather in the gauging procedure itself. Analogous 
subtraction terms arise in the three-nucleon problem whose strong interaction equations have no subtraction terms 
{jJU- Similar subtraction terms will arise in the gauging of -Fi and F\. The graphical representation of A4 M and V t M 
is given in Fig. |3|. 

Eq. ( [f0| ) can be used to determine all the possible electromagnetic transition currents of the irNN system. An 
especially interesting use of Eq. ([lo|) is to study the electromagnetic properties of bound states of the ttNN system. 
It is certainly expected that the strong interaction irNN model under discussion admits a bound state corresponding 
to the physical deuteron. In this case a solution will exist to the homogeneous version of Eq. (^): 

$ = Vg t <f> (42) 



G 



where $ = ( $n <£>a &d) T ■ Here &n is the usual deuteron vertex function describing the d — * NN transition, 
while $a and <E>d are somewhat unusual in that they describe transitions to clustered ttNN states: d — > (irN)N and 
d — > (NN)tt respectively. Comparing Eq. (p7j) and Eq. (p?) it is seen that Q has a pole at the deuteron mass M4: 



as P^ Mj (43) 



where P is the total four-momentum of the system, and where ^ satisfies the equation 

* = GaVt^. (44) 

Clearly is related to $ by the equation \P = ^ ^ <f> so that either of the equations (fl2|) or ( fii]) can be used 

to determine 'J. Taking the left and right residues at the deuteron bound state poles of Eq. ( f40| ) we obtain the bound 
state electromagnetic current 

f = * (M" + Vf) * (45) 
which describes the electromagnetic properties of the deuteron whose internal structure is described by the present 



ttNN model. Eq. ( 45 ) provides a very rich description of the internal electromagnetic structure of the deuteron with 
all possible meson exchange currents being taken into account in a gauge invariant way. In view of the accuracy of 
this model which is based on meson and baryon degrees of freedom, a comparison of the deuteron electromagnetic 
form factors (easily extracted from j M ) with experiment should prove to be very interesting. 
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